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Abstract 

In the framework of the Ginzburg-Landau harmonic potential approxima- 
tion, we present a possible modeling of the time-dependence of the frequency 
of the order parameter mode suitable to account for the formation of corre- 
lated domains in non-instantaneous symmetry breaking phase transitions. An 
interesting spectrum of possibilities for the size and the life-time of these do- 
mains emerges, which appears to be consistent with the conclusions of other, 
differently grounded analysis. 
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I. INTRODUCTION 

In some recent papersEM much attention has been devoted to phase transitions in the 
so-called range of criticality, namely in the temperature range spanning between the critical 
temperature Tc, at which the phase transition "starts", and the Ginzburg temperature Tq, 
at which the broken symmetry phase becomes maximally stable. Long string formation in 
such a region has been hypothesized and its consequences have been investigatedB. Beyond 
the detailed phenomenological picture and the numerical simulations, the dynamic formation 
of non-homogeneous structures in the critical region seems to require further study. Our aim 
in this paper is to describe indeed a possible mechanism for the formation of domains in the 
process of non-instantaneous phase transitions characterized by time dependent order pa- 
rameter. By resorting to known results in the quantum field theory (QFT) with spontaneous 
breakdown of symmetry (SBS), we also discuss, on the basis of a microscopic analysis, the 
assumption of temperature and time dependence of the phenomenological Ginzburg-Landau 
(GL) chemical potential. 

To this aim let us first shortly summarize some known features of the spontaneously 
broken symmetry theory with non-constant order parameter. 

We are interested in the case where a specific phase of the system at a certain time is 
characterized by a specific value of the order parameter. Changes in the order parameter, 
considered as a function of the temperature or of the time, or of both of them, are thus 
associated to a change of the system phase, i.e. to the process of phase transition. Here 
we closely follow reLu. The microscopic analysis may start considering, for simplicity, a 
quantum real scalar field (pix) at T 7^ with Lagrangian and Hamiltonian given by 

£ = ^9^09^0 - ^/iV - ^A0\ A>0, (1.1) 



H = J d^'-'xllidocPr + \{dcj>Y + i/iV' + \^'] , (1.2) 

respectively, where D is the space-time dimensionality and x = (x, t). The equation of 
motion is 

i-d'-fi')<P = \<P' , (1.3) 

where the conventions g^^ = —g^^ = 1 and c = h = 1 have been adopted. In the variational 
approach, the free energy functional is given by the Bogoliubov inequality! 

F<Fi = Fo+<H-Ho>o , (1.4) 

where F = —KT lnTr[exp(— /5if)] , with [3 = is the free energy and the symbol < >o 
denotes the statistical average: < A >o= {^)TT[exp{—(3Ho)A], with Z = [exp{—(3Ho)]. Hq 
is the trial Hamiltonian. One introduces the variational parameter /xq which in general can 
depend upon x and P and acts as the mass for the quantum field p. For a discussion on the 
relation between p and the field see ref.i where the detail of the variational computation 
is presented. Here we only recall that we have (p = p + v, with the c-number v being the 
non-vanishing expectation value (the order parameter) for the (p field: 
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< O|0(x)|O >= V 



(1.5) 



Eq. ( |1.5|) implies /x^ < and expresses the condition of spontaneous breakdown of symmetry. 
A specific physical phase of the system corresponds to a specific value of v. In the present 
paper we are interested in the case of a non-constant order parameter: v = v{x,l3). In such 
a case, Eq.(|T]3) tells us that the ground states or vacua may be labeled by the values of the 



order parameter: |0 >= |0 >^[-rcl. Thus different vacua (i.e. different phases) correspond 
to different values of the order parameter: |0 >„7^ |0 >^', with v' = v{P'), v{P) ^ f (/?') 
and /3 7^ In order to study the process of phase transition one thus searches for the 
equation describing the variations of the order parameter f(a;, /3). Such an equation may 
be understood as describing the "motion" over trajectories in the "space of the vacua" (the 
phases) of the system (also called the space of the representations of the CCR). Let us then 
observe that eq. ( |1.4| ) gives (we are assuming real f(a;, /?)) 



Fi(^;(x,/5)) = y ip-^x\{d^vf - L,^^{v{x,(5),d^v[x,m (1.6) 

with 

l^efJ = \d'vd,V-V,ii{v) , (1.7) 

UesA^) = Ucliv) + FM) -h<p'>l , (1.8) 

U,l{v) = \^^\' + \\v' . (1.9) 

Note that Ueff{v) includes not only the classical terms but also the quantum and the 
thermal contributions. We remark that the free energy Fi plays the role of the energy 
for the c-number field v{x,(3) whose Lagrangian is Ceff'- in this way the dynamics for 
the quantum field manifests as the dynamics for the classical field describing the order 



parameter f(x,/3). Eq. ( p..6|) and Ueff{v) are the generalized Ginzburg-Landau functional 
and potential, respectively. The Euler-Lagrange equation of motion for the field v{x,P) is 
deduced from C^ff through the action principle: 

{-d^ + m^)v{x, (3) = Xv%x, (3) , (1.10) 

with 

= l/i^l -3A < >o (1.11) 

and therefore m? is temperature dependent: rn^ = m?{(3). The variational equation |^ = 
gives the variational mass parameter 

/ig(x, /?) = + 3Af 2(x, /3) + 3A < p2 >o . (1.12) 



Eq.( |1.10D describes the evolution of the order parameter, namely the evolution of the sys- 



tem in the "space of the representations of the CCR", each representation being labeled 
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by v{x,P), and therefore the evolution through the different phases of the system (phase 
transitions). Eg . ( p. . is a consistency relation between the "inner" dynamics described 
by £ (eq.(|l.l|)) and the dynamics in the space of the representations described by £e// 
(eq.(|3))- Notice that the " mass term" in eg . (11.101) has the "wrong" sign with respect to 
a Klein-Gordon eguation for a physical particle field: — is interpreted as the GL chem- 
ical potential and negative chemical potential signals the realization of SBS. Temperature 
dependence in the mass term is usually introduced by hand. Here, through eg. ( |1 . 1 1| ) we see 
the microscopic origin of the temperature dependence of the chemical potential in the GL 
potential. 

In the case of one space dimension one may obtain the stationary kink solution!. When 
a gauge field is present and is a complex field one has the vortex solution in two space 
dimensions. In theories with S0(3) and SU(2) symmetry groups and three space dimensions 
one has the monopole and the sphaleron solutionJl^. A popular framework where the for- 
mation of defects in the process of non-eguilibrium symmetry breaking phase transition is 
studied is the Kibble-Zurek one0. 



II. THE GINZBURG REGIME 

Having presented the essential outline of the framework of our discussion, we now turn 
our attention to the transition processes occurring in a finite span of time (see, for exampleSi) 
in which new phenomena, such as the formation of exponentially decaying strings, are ex- 
pected to occur. In these processes the transition starts at the critical temperature Tq and, 
after a certain lapse of time, the maximally stable configuration is attained at the so-called 
"Ginzburg temperature"^ Tq {Tq < Tc). Between Tc and Tq the system is said to be in 
the critical regimeS. 

Our purpose is to model the time dependence of m?{(3) during the critical regime evo- 
lution, i.e. for transitions lasting a finite time interval where the formation of domains is 
allowed. 

Following ref.0, we will focus our attention on the harmonic limit which is obtained 
by considering only the (linear) l.h.s. of eg.( |1.10|) . Of course, this introduces a strong 
constraint and it does not accurately describe the behavior of the system. However, it gives 
enough reliable information on the critical regime behavioJll. We notice that in such an 
approximation /Iq = —m?. 

Let us now consider the expansion of the w-field into partial waves 

v{x, (3) = Y: {Mt, /5)e^''-^ + ulit, /3)e-^'^-^} . (2.1) 



Then, in the harmonic potential approximation, the GL eguation (|1.10| ) gives for each k- 
mode {k = -\/k^): 



Uk (t) + (k^ - m^)uk = 0, (2.2) 

i.e. it leads to the eguations for the parametric oscillator modes (see also0). Let us 
denote the freguency of the fc-mode by 



Mk{t) = sjk^-m^{f3,t). (2.3) 



4 



which is required to be real for each k (in this paper we are not interested in fact in unstable 
modes corresponding to imaginary frequency; for a recent discussion on these modes and 
spinodal decomposition see, e.g., ref.i and references therein quoted). In full generality, we 
assume that as well as f3 may depend on time. We remark that the reality condition on 
Mk{t) turns out to be a condition on the fc-modes propagation. The reality condition in fact 
is satisfied provided at each t, during the critical regime time interval, it is 

k2>m2(/3,t), (2.4) 

for each fc-mode. Here and in the following, for shortness we use the notation j3 = (3{t). Let 
t = and t = T denote the times at which the critical regime starts and ends, respectively. 

For a given k, eq.(|2.4|) will not hold for any t: of course, it holds up to a time t = Tk after 
which is larger than fc^. Such a represents the maximal propagation time of that k 
mode. The value of is given when the explicit form of is assigned. In the following we 
consider two possibilities for the time-dependence of m(/?, t). The first one can lead to large 
correlation domains, and thus possibly excludes defect formation. The second one does not 
allow very large domains and therefore it allows defect formation. Our first model choice is 

a) m2(/?,t)=m2(e^^(^«'*)-l), (2.5) 

with t = assumed to correspond to the minimum of (we assume /(O) = in the 
example of Figure 1 and in the following). Notice that although at the transition temperature 
infinite correlation length is allowed, the corresponding mode has only a limited time for 
propagating. So the formation of domains, i.e. the 'effective causal horizon can be inside 
the system (domain formation) or outside (single domain) according to the time occurring 
for reaching the boundaries of the system is longer or shorter than the allowed propagation 
time. So, the dimensions to which the domains can expand depend on the rate between the 
speed at which the correlation can propagate, at a certain time, and the correlation length 
(Afc oc (m(/5,t))~^) at that time. 

In the present case, each fc-mode can propagate for a span of time < t < r^. From 
eqs.( p.4|) and ( p.5|) we obtain: 



f{Pin),n) = In ( '^'^^""' j oc In (I) , (2.6) 

where and Sq are the fc-mode energy for non-zero and zero k, respectively. The equilibrium 
time at which T = To is t > Tk for any k. 
A second possibility to model is: 

b) m2(/3,t) = m2e2^(^W'*) . (2.7) 

In this case, a cut-off exists for the correlation length, L oc mg ^ (in the example of Figure 2, 
it corresponds to the time t = 0). In this case the propagation time is implicitly given by: 



f{P{n),n) = ln{ — \ (xln - , (2.8) 
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f{f3{Tk), Tk) resembling the commonly called string tensionEl. In eq. ( |2.8|) ^ is the correlation 
length corresponding to the k-mode propagation. In this case the reality condition acts as an 
intrinsic infrared cut-off since small k values are excluded, due to eq.( p.4| ), by the non-zero 
minimum value of m^. This means that infinitely long wave-lengths are actually precluded 
i.e. only domains of finite size can be obtained. Finally (phase) transitions through different 
vacuum states (which would be unitarily inequivalent vacua in the infinite volume limit) at a 
given t are possible. This is consistent with the fact that the system is indeed in the middle 
of a phase transition process (it is in the critical regime )lli. At the end of the critical regime 
the correlation may extend over domains of linear size of the order of oc {m{f3o , t))~^ . 

In conclusion we see that the model choice a) differs from the model choice b) in the 
fact that the case a) allows the formation of large correlation domains (infrared modes are 
allowed); in the case b), on the contrary, only finite-size domains can be formed. Due of the 
fact that defects may be thought as the normal symmetric phase region trapped in between 
correlated (ordered) domains, we see that model a) may exclude defect formation if, as 
mentioned above, k-mode life-times are longer than the ones needed to reach the system 
boundaries. 

We may then specify the function / by assuming that it is positive, monotonically growing 
in time from t = to t = r, at which the Ginzburg temperature is reached. 
A possible analytic expression for f{P{t),t) is 

where a, b, c are (positive) parameters chosen so to guarantee the correct dimensions and 
the correct behavior in time and in temperature (see below for the relations between the 
various time scales and the ratios among a, b,c): T decreases from Tc to Tq as time grows 
from t = to t = t; f{j3{t),t) is positive in the critical region (and even for t > t, i.e. in the 
SBS region). The equilibrium time scale is given by = |. Notice that in order to obtain 
this result, the variations in time of (5{t) have been assumed small, according to the picture 
of a slow (non-instantaneous) transition (here we refer to slow transitions as those for which 
r is large; for fixed parameter a, large r also means large tq, see below). 

The behavior of m?{(3,t) corresponding to ( |2.9| ) is shown in Figure 1 (case a) and in 
Figure 2 (case b). The maximum of m?{(3) has to be identified with the minimum of the 
potential = = — m(/9)^v(a;, /3) + Af ^(x, /3) (cf. eq. (|l.ll|) ), so it corresponds to T = Tq- 
The negative time region in the Figures 1 and 2 corresponds to T > Tc, so there / should 
be negative. 

Finally, the present treatment allows to recover the known results on the number of 
defects, i.e. n^ef oc {jq/tq)^ (see ref.i, for a recent review). In fact, recalling that the 

equilibrium time r is given by r = \Jc/b it is possible to introduce the time scales: tq = c/a\ 
and To = aA/6, with A arbitrary constant and a, 6, c introduced in eq. (pl9|) . Thus r = ^tqTq, 
which also says that for fixed a (and b) large r also means large tq. We now observe that, 
at the first order approximation, one has 

e^(-) - g{x) ^ e^'^^-") - g{xo). if = 0, (2.10) 

so, in our case, e^-^*^"^^ = 2/(r) + e^f^'^^ 
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The number of defects can be finally obtained. Our choice is to consider the defect 
number at Tq, namely when the system gets enough stabilized (in some sense we count the 
defects that are still present at the end of the critical regime). In the case a) it is: 

Udef oc m^{l3G,r) 2mlf{T) = ml t/Xtq oc \IJt^/tq). (2.11) 

Similarly, in the case h): 

ridef oc rn?{(3G, t) ~ mQ(2/(r) + 1) ^ ttIq t/Xtq for large r (slow transitions) (2.12) 

As observed above, for fixed a (and b) large r also means large tq. We note that the form of 
is a good approximation at the Ginzburg temperature in our models, since, as already 
observed, at Tq it corresponds to the minimum of the potential, which is where we want 
to know the defect number. It is in fact interesting that Zurek's results^! are there also 
recovered. 



III. FINAL REMARKS 

It is interesting to write M|(t) as 

M,2(A,(t)) = M,2(0) e-^^^W, (3.1) 

sinh(/(r,)-/(/3(t),t))^ ^22) 
smh/(rfc) 

where /(O) = has been used. Eq.(g]2|) shows that Afc(t) > for < t < , Afe(O) = 
and Afc(rfc) = oo. Since Mfc(Afc(rfc)) = we see that Ak(t) acts as a life-time, say with 
Afc(t) oc Sk, for the fc-mode: each fc-mode "lives" with a proper time Sk, i.e. it is born when 
Sk is zero and it dies for Sk — > oo. The "lives" of the fc-modes for growing k are shown in 
Figure 3, where Ak(t) is plotted versus t for different values of k. The blowing up values are 
reached in correspondence of the r^, values for each fc-mode. Only the modes satisfying the 
reality condition are present at a certain time t, being the other ones decayed. In this way 
the causal horizon sets up. Figure 3 and Eqs. ( p.lD and ( p.2| ) show that modes with larger 
k have a longer life with reference to time t. 

Since, on the other hand, longer wave-lengths correspond to lower k, we see that domains 
with a specific spectrum of /c- modes components may coexist, some of them disappearing 
before, some other ones persisting longer in dependence of the number in the spectrum of 
the smaller or larger k components, respectively. In general, the boundaries of larger size 
domains are thus expected to be less persistent than those of smaller size domains. This fits 
with the observations that " critical regime has little effect over the small scale dynamics" , 
thus allowing the survival of localized defects (such as vortex strings). 

In conclusion, we have studied, in the framework of the GL harmonic potential approx- 
imation, the non-equilibrium symmetry breaking phase transition processes occurring in a 
finite time interval. We have presented a possible modeling of the time-dependence of the 
frequency of the order parameter modes suitable to account for the formation of correlated 
domains. An interesting spectrum of possibilities for the size and the life-time of these 
domains emerges, which appears to be consistent with the conclusions of other, differently 
grounded analysis. 
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FIG. 2. Case b) Behavior in time of ■m?{(3,t) 
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FIG. 3. Life-time for different values of k 
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